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Abstract 

We construct crude estimates for non-optimality of quantum measure- 
ments in terms of their violation of Holevo's simplified minimum-error 
optimality conditions. As an application, we show that a modification of 
Barnett and Croke's proof of the optimality conditions yields a convergent 
iterative scheme for computing optimal measurements. 
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1 Introduction 



The minimum-error quantum detection problem arose in the 1960's in the design 
of optical detectors [1] and has been of recent importance in the subjects of quan- 
tum information [31 dj [5j and quantum computation [51 [71 [HJ [SI [13 fTT] : 

If an unknown state pu is randomly chosen from a known ensemble 
of quantum states, what is the chance that the value of k will be 
discovered by an optimal measurement? 

Barnett and Croke Jjj have recently provided a simple operator-theoretic 
proof of the necessity of the standard Yuen-Kennedy-Lax & Holevo (YKLH) 
optimality conditions |131 114j for the minimum-error quantum detection prob- 
lem. Their proof may be shortened, since Holevo [15] had previously shown that 
an intermediate step of their proof (positivity of the operators Gj defined by 
equation (10) of |12| ) provides a simplified necessary and sufficient condition for 
minimum-error quantum detection. 

1.1 Results 

This note gives a more robust version of Holevo's simplified optimality condition 
(condition II of Theorem [21 below), by estimating non-optimality in terms of 
quantitative violation of this condition. These bounds are used to show that the 
perturbative method of Barnett and Croke may be converted into a convergent 
iterative algorithm for computing optimal measurements, adding to the list 
[El [ITl [iHl [T9j of algorithms for this purpose. This iteration converges even 
for countably-infinite ensembles in an infinite-dimensional Hilbert space. 

2 Conditions for minimum-error quantum dis- 
crimination 

A precise description of the minimum-error quantum measurement problem is 
given by: 

Definition 1 Let 8 = {pk\keK ^'^ ensemble of mixed quantum states Pk, 
which are represented as positive semidefinite operators on a Hilbert space Ti 
normalized by a-priori probability: Tr pk = Pk with ^^Pfe = 1 • The support 
supp {£) is the closure of the span of the ranges of the pk- A positive operator- 
valued measurement (POVM) is a collection of positive semidefinite oper- 
ators {Mk} satisfying ^^M^ = 1. The corresponding Lagrange operator is 
given by 




(1) 
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The minimum-error quantum discrimination problem JT] consists of 
finding a POVM maximizing the success probability 

Psucc {{Mk}) - Tr J^MkPk = TtL (2) 

k 

of correctly distinguishing an element blindly drawn from the ensemble £ . (We 
will often abuse notation by writing Psucc {Mk) instead of P^ucc {{Mk})-) 

Holevo's simplified optimality conditions are given by property II o|3 

Theorem 2 (Holevo [TS], Yuen-Kennedy-Lax [13j, ) Let{Mk}f.^^ ^ be 
a POVM for distinguishing the ensemble £ . Then the following are equivalent: 

L {Mk} maximizes 

Psucc ■ 

IL [l + L'^ /2 > pk for all 

III. There exists a self-adjoint operator G satisfying G > pk and {G — pk) Mk = 
for all k. 

Furthermore, under these equivalent conditions L = = G, and L is the unique 
self-adjoint operator of minimal trace satisfying L > pk for all k. 

The above optimality conditions were first proved in the infinite-dimensional 
case by Holevo, since earlier proofs worked only in finite dimensions. The in- 
equalities in properties II-III use the standard order on self-adjoint matrices: 
A> B iS A — B is positive semidefinite. The LHS of condition II is commonly 
referred to as the real part: 

Re (L) := {L + L^) /2. (3) 



3 Mathematical background 

Definition 3 Let A be a self-adjoint operator on a Hilbert space Ti. with spectral 
decomposition A = ^^XkTlk- The positive part of A is given by 

[A]^=Y.XkIlk. (4) 

Afc>0 

The positive projection is given by 

X+{A)^Y.^k. (5) 

Afc>0 

The trace norm of an operator B -.Ti ^ Ti is given by \\B\\^ = TrVB^B. The 
operator norm is given by \\B\\^ — sup||^||^]^ II^V"!!- 

^Another interesting optimality condition is given by Theorem 3 on page 39 of 12UI . 

^Earlier formulations of condition II |13l I14| were that L = L' and L > p^. for all k, 
equivalently stated as L = and Re (L) > p^.. (The self-adjointess condition is redundant 
in the latter form.) 
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We collect some simple mathematical facts. We will frequently use the in- 
equalities 

|TrA| < II All, (6) 
\\BC\\,<\\B\\,\\C\\^, (7) 

which may be found in [21] . For positive semidefinite operators Pi , P2 > such 
that P1P2 is trace class, one has 

TrPiP2>0, (8) 
with equality iff P1P2 = T3] and 

Ai>A2^ CUiC > CU2C (9) 
for all operators C and self-adjoint Ax.Ai- 

4 Estimates of near- and non-optimality 

Our next goal is to strengthen condition II of Theorem [2] by giving quantitative 
bounds in the case that condition II fails to hold. As a first step, note that in 
the finite-dimensional case if 

Re {V) >pk-a (10) 
for some scalar a > 0, then by inequality ^ 

Psncc [Mk ) - Tr Re (L) = Tr ^ Re (L) M°p* (11) 

k 

> Tr ^ [pk - a) M°P* = P,„,e (mD ^ « dim (12) 

k 

where M^''* is some optimal POVM. 

In order to control dimensional factors (and to consider ensembles on infinite- 
dimensional Hilbert spaces) it is useful to introduce the following concept: 

Definition 4 Let £ — {pk} be the ensemble of definition]^ and let p G [0, 1]. 
The p-dimension dinip {£) is the minimum dimension of a subspace A for 
which 

Y^Wii -nA)pkh<P, (13) 

k 

where Il\ is the orthogonal projection onto A. 
Remark: Note that the inequality implies that 

Tr(i-nA)^Pfc < ||^(i-nA)pfe||^ <^||(i-nA)pfc||, <p. 
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Lemma 5 For fixed £ , the function p t-^ dinip (£) is finite on (0, 1] and monotonically- 
decreasing on [0, 1]. 

Proof. The monotonicity of p '—^ dinip (£) is immediate from the definition. To 
prove finiteness for p > 0, take spectral decompositions pk = ^^A^^ IV'm) {ipke\- 
For a finite subset S of the (fc, £), let lis be the projection onto the linear span 
of the ipke with {k,i) G S. By the triangle inequality 

^||(i-n5)pfc||i <5]||(i-ns)Afc,|V^fe£)(V'w||li < J2 

k ki (k,l)<^S 

Since X^i ~ 1, we may take a finite subset S of the (fc,^) such that the 

right-hand side may be made smaller than p. ■ 

We may now state a robust version of Theorem [2l 

Theorem 6 Let {Mk} be a POVM for distinguishing £, let L — ''^^M^pk, and 
let {MjJ?^*} be an optimal measurement. Then 

1. Assume that a > is a scalar such that 

Re {L) >pk~a (14) 
for all fc. Then for p e [0, 1/4) 

Psucc (Mk) > Psucc (Afro - " dimp {£) - Ap. (15) 

2. Suppose that Re (L) ^ pi for some £. Then 

Psucc (Mk) < Psucc (M°^") - (Tr {[pe - Re (L)]+))' , (16) 
where [•]_^_ is the positive part, defined in definition\^ 

4.1 Discussion of Theorem [6] 

The small-a case of Part 1 addresses the case where {Mk} nearly-satisfies con- 
dition II. In particular, (fT^ implies that Psucc (Mk) > Psucc {M'jl'^^) — e if 

£ — ip 

a < sup (17) 

pe[o,£/4] climp (t) 

The following example shows that the dependence of this expression on £ may 
not be removed except (in the finite-dimensional case) by introducing dimen- 
sional factors: 

Example 7 Let m be a positive integer, and let £ be the m- state ensemble on 
C™ defined by pk = \k) (fc| /m. Set Mk = |fc -I- 1) (fc -I- 1|, using addition mod m. 
Then one has Psucc (Mk) = and Psucc (^f^*) = 1, but inequality ( |i^[ ) holds 
for a — 1/m, which approaches as m — > oo. 
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4.2 Proof of part 1 of Theorem [6] 

Proof. Let n be an orthogonal projection, and set = 1 — 11. Then 

-Psucc (Mfc) = Tr (n Re (L) H) + Tr (n^ Re (L) U^) . (18) 
Using equations ©-([H]) to estimate the first term, 

Tr n Re (L) H = Tr ^ Re (L) x nA/°P*n 

fe 

>Tr^(pfe-a) xnM°P*n 

= Psucc (M°P') - a Tr (H) + Tr {Up^U - p^) M°p') 

> Psucc (M°P') - a Tr (n) - ^ \\pk - npkU\\^ . (19) 

But 

J2\\pk- ^PkU\\, = J2 W^^Pk + PfcH^ + n^PfeH^ II 1 

k k 

<zY.\\T\^Pk\\^. (20) 
Using ([7|) to estimate the second term of ([TSl) . 

|Tr(n^Re(L)n^)| < i ||^n^pfeii4n^ + n^MfePfeH^ ||^ 

<EI|n^Pfe|li (21) 

Putting (UHl) - (HI]) together gives 

Psucc (Mfe) > Psucc (Af°P*) - a Tr (H) - 4^ ||n^pfe || ^ . (22) 

The bound ([15]) follows by picking 11 to minimize Tr (11) when the last term of 
(|22)) is constrained to be less than p. (By Lemma [5] such 11 of finite rank always 
exist.) ■ 

4.3 Proof of part 2 of Theorem [6] 

Definition 8 Let {Mk\ he a POVM for distinguishing the ensemble £ of defi- 
nitionl^ let X < 2 x 1 be a positive semidefinite operator on Ti, and let i G K . 
Then the Barnett -Croke modification of {Mk} is defined by 

Mk {X, ^) = (1 - X) Mk {1-X)+ Ski (2X - X^) . 

Remark: Note that since < 2X - for < X < 2 x 1, for each £ the 
set {Mk{X,£)} forms a POVM. Barnett and Croke [12] considered the case 
X = e lip) ('01, where ip is a unit vector satisfying the eigenvalue equation 

{pe-Re{L))\yj)=-X\ij), 
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with A > 0. They showed that 

Psucc (Mfc(X,^)) = 2A>0. 



d 
de 



e=0 



In order to complete the proof of part 2 of Theorem [6l it suffices to tm'ii this 
perturbative argument into an estimate. 

Proof of part 2 of Theorem [6l Let 11+ be the positive projection ([5]) 

n+ = x+(Pf-Re(L)). (23) 

Then for a e [0,2], 

Psucc (Mfe {all+,e)) = F.ucc (Mfc) + 2a Tr [{pe - Re (L)) x n+] 

- Tr (U+pe) + Tr ^H+A^n+pfe (24) 
> Psucc (A4) + 2a Tr ([p, - Re (L)]+) - (25) 

where we have used cychcity of the trace and® — ©■ 
Note that if Tr {[pi - Re (L)]^) > 1 then 

P,ucc (Mfe {n+,e)) = Fsucc (Mfe) + 2 Tr ([p, - Re (L)]+) - 1 > 1, 
giving a contradiction. In particular, we may set 

a-Tr([p,~Re(L)]^) e [0,1], (26) 
maximizing the RHS of ([25]) over a E [0, 1]. This gives 

P.UCC {Mk (an+, £)) > P,„ec (A4) + (Tr ([p, - Re (i)]+))' . (27) 



5 Barnett-Croke iteration 

In this section we show how to convert Barnett and Croke's perturbative proof 
into an algorithm for computing optimal measurements. Although the success 
rate of poorly-chosen iterations might fail to actually converge to that of an 
optimal measurement H the following sequence does not exhibit this malady: 

Definition 9 Let {M^} be a POVM for distinguishing the ensemble £ of defi- 
nition]^ and chose t to maximize 

a = Tr [p, - Re (L)] (28) 



•^In is asserted in |17| that the algorithm of ,16. suffers this fate. 
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Then the iterate of {Mk} is the POVM 

M+ = Mk {ax+ [pi - Re (L)) , I) , (29) 

where [•J^ and x+ are defined m Q — (0). For a given measurement {A/^"-*}, 
recursively define the iterative series {Affc"^}r!>i 

Mf."+^) - (a^"^)^. (30) 

Remark: An index I maximizing (|28p exists using minimax principle (Theorem 
XIII.l of [22) and the fact that Tr^p£ = 1. 

The proof of part II of Theorem [6] actually proved the following stronger 
result: 

Theorem 10 The above iteration monotonically increases success rate. In par- 
ticular, for an arbitrary POVM {Mk} the set |A/^} is a well-defined POVM, 
and 

Psucc (A/+) > Psucc {Mk) + max (Tr {[pt - Re (i)]+))' . (31) 
We now show that the iterative scheme of definition[9]approaches optimality: 

Theorem 11 Let M'j^^ be an arbitrary starting POVM for the iterative series 
((23). Then 

hm PsuUmI^^) = PsucrXM°J'), (32) 

n— *oo 

where Af^^* is an optimal measurement. 

Proof. Let e > be arbitrary. We seek an > such that 

n>N^ Psucc (Af^"^) > Psucc (Af°P') - £. (33) 



Set 



L(")=^A^(^ 



in) 



Pk- 



By equation (jl7p and the monotonicity of n ^ Psucc (m^"^), it sufHces to find 
& n < N such that 

Re >Pe-A (34) 

for all £, where A is any real number satisfying 

0<A< sup ^ (35) 



^Faster convergence can be obtained by replacing a by /3 in equation I I29I I , whore /3 € [0, 2] 
is chosen to maximize Psucc (^^k^' which is quadratic in /3. 

^In finite dimensions, one may take A = e/ dimH. < e/ dim (supp (£)), corresponding to 
p = 0. 
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We claim that N = A ^ suffices. 
Assume that 



maxTr 

e 

for all n < A^. By Theorem [TUl 



Pi - Re 



> A 



^succ ( ^'^k 



> TV X > 1, 



yielding a contraction. 
It follows that 



maxTr 

e 



< A 



Pi - Re 

for some n < N. The inequality ([M)) follows from the observation that 

A < Tr {[A]^) X 1, 

for A = p£-Re(L(")). ■ 



6 Conclusion 

Using non-optimality estimates in terms of quantitative violation of Holevo's 
simplified optimal measurement condition, we have converted Barnett and Croke's 
perturbative proof into a conceptually-simple iterative scheme for computing op- 
timal measurements. This iteration approaches the optimal success rate even 
in the case of infinite-dimensions and infinite ensemble cardinality. It would be 
interesting to try to improve the non-optimality bounds of Theorem |6l and to 
study the convergence rate of this iteration in more detail. 
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